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ABSTRACT 

We present a high spatial resolution optical and infrared study of the circumnuclear region in Arp 
220, a late-stage galaxy merger. Narrowband imaging using HST/WFC3 has resolved the previously 
observed peak in Hq;-|-[NII] emission into a bubble-shaped feature. This feature measures 1.6” in 
diameter, or 600 pc, and is only 1” northwest of the western nucleus. The bubble is aligned with the 
western nucleus and the large-scale outflow axis seen in X-rays. We explore several possibilities for the 
bubble origin, including a jet or outflow from a hidden active galactic nucleus (AGN), outflows from 
high levels of star formation within the few hundred pc nuclear gas disk, or an ultraluminous X-ray 
source. An obscured AGN or high levels of star formation within the inner ~100 pc of the nuclei are 
favored based on the alignment of the bubble and energetics arguments. 


1. INTRODUCTION 

Galaxy mergers are integral to understanding galaxy 
evolution. Simulations suggest that major mergers cause 
large-scale tidal gas inflows, which may fuel nuclear 
starbursts and active galactic nucleus (AGN) activity 
(Barnes & Hernquist 1996 |Mihos & Hernquist 1996). 
In later stages, mergers may host young, obscured AGJN, 
which, simulations predict, are gradually uncovered as 
AGN feedback and high levels of ci rcumnuclear star for¬ 
mation expel the surrounding g as ( Barnes fc Hernquist 
19961 IMihos & : Hernquist||1996[ |Di Matteo et al .||2L)U51 

IHopk: 


Hopkins et al. 2005 IHopkins & Hernquist||2006p. This 


nuclea r activity can, in turn, fuel large-scale galactic out¬ 


flows (iHeckman 2003 Veilleux et al. 200^, which have 
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Large-scale outflows in the outer parts of galaxies have 
been extensively studied. However, the advent of high 
spatial resolution observations has prompted more recent 


cape of these outflows (Rupke & Veilleuxl2013a 

Emonts 

et al. 

2014 

Sakamoto et al.[ 

2014). When observed 


with suthcient spatial resolution, an outflow still within 
a galaxy often has complex superbubble or filamentary 
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gion interacts with the cooler interstellar medium (ISM). 
However, these well-shaped structures have thus far not 
been observed in the turbulent, clumpy ISM of a late- 
stage merger. In addition, while both the timescale and 
impact of AGN feedbac k into the host system have been 


calculated theoretically (Hopkins et al.|2005 Hopkins & 
Elvis|[2M0|, both are still observationally uncertain. In- 


vestigatmg an individual late-stage merging system that 
exhibits both an outflow and a possible obscured AGN 
may provide clues to this process. 

Arp 220 is the closes t ultraluminous infryed galaxy 
(ULI RG, Lju >10^^ Lq, [S anders et al.|l988 Soifer et al. 


1987). The system has long been recognized as a late- 


stage merger, a s evidenced by extended tidal ta ils visible 
in the optical (Arp 1966 Surace et al. 2000), infrare d 
( Joseph fc Wrightf .985), and in HI ( Hibbard et al.|2000 ), 
as shown in figure It also hosts a large-scale outflow 
that extends far from the galaxy. Arp 220 ha s a molecu¬ 
lar gas mass of ~10^° Mq (Scoville et al. 1997) and a total 


star formation rate, calculated using the far-infra red lu¬ 
minosity, of 240±30 Mq yr“^ ( Farrah et al.|200'3 |. Most 
of the molecular gas and star formation are concentrated 
within a central few hundred pc nuclear gas disk, which 
also hosts a double nucleus. A possible AGN has been 
suggested (Paggi et al. 2013) but is still debated. At a 
distance of 77 iVipc, its close proximity offers a unique 
laboratory to investigate the interaction of its outflow 
with the ISM at small spatial scales, and to determine 
how energy propagates out of the interior of such a sys¬ 
tem. A more detailed background on the structure of 
A rp 220 is present e d belo w. 


Heckman et al. (1987) observed large-scale (10 kpc) 
bipolar lobes extending from the galaxy via narrow- 
band Ha imaging, which they interpret as a superwind. 
Further spectroscopic study reve aled these lobes to b e 
superwind-driven bipolar bubbles (Heckman et al.|1990). 
This superwind is also observed in soft X-ray einission 
(McDowell et al.||2003). 

In addition to the large-scale wind, smaller-scale gas 
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Figure 1. Left: HST ACS image of Arp 220 as first presented in [Wilson et al.| (|2006|l, rotated to adjust for the 77° position angle. The 
blue image is taken with the F435W filter and the red with the F814W niter (green in a linear combination of the red and blue filters). The 
gray box indicates the FOV of the zoomed right panel. Right: Our new data taken with HST WFC3. Blue is the F467M blue continuum 
image, green is [OIII], or the FQ508N image with the blue continuum subtracted, and red is Hq:+[NII], or the F665N image with the red 
continuum subtracted. North is up, east is left. 


dynamics have b een investigated in more recent opti¬ 


cal IFU studies. Colina et al. (2004) found large-scale 
plumes and lobes, up to 8 kpc in length and aligned with 
the superwind, which are consistent wit h high-velocity 
shocks expanding in a neutral medium. [Arribas et al. 


(2001) studied the inner region of the system within 2 


kpc of the nuclei and found three gas components. One 
component is coupled to the central molecular gas disk. 
The other two components are associated with a biconal 
outflow or superwind. 

This outflow, whether star formation- or AGN-driven, 
almost certainly originates within the nuclear region 
of the galaxy. At the center of Arp 220, two dust- 


ble in 

the IR ([Scoville et al. 

1998), sub-mm ([Scoville et al. 

1997 

Downes & Solomon 

19981 Sakamoto et al.l 1999 

Scovil 

e et al.||2014|), and radio (Baan Haschick|1995D. 

Sakamoto et al. {1999 

) resolved the large circumnuclear 


cleus is embedded in its own gas disk, each roughly 100 
pc in radius, and these disks are counterrotating with 
respect to each other and inclined such that their near 
sides are to the south. The position angle (P-A.) of each 
of these small disks is 52° and 263° for the eastern and 
western nuclei, respecti vely, while that of the larger cir- 
cumnuclear disk is 25°. Scoville et al. (2014) calculated 
the interstellar medium (ISiVlj masses of each nuclear 
gas disk, finding a mass of 1.9 (4.2) xlO® Mq around 
the eastern (western) nucleus, concentrated within a ra¬ 
dius 69 (65) pc of the nucleus. These small gas disks 
are surrounded by a larger gas disk ^1 kpc in radius, 
whic h rotates around the d ynamical center of th e sys¬ 
tem (iSakamoto et al. 1 19991) . iMundell et al.l (|2001|) used 


observations of HI in the circumnuclear gas disks and 
comparison with simulations to estimate that the merger 
started about 700 Myr ago, with the second, most recent 
starburst starting in the western nuclear gas disk more 
recently than 10-100 Myr ago. 

Within the inner 50 pc of each nucleus, the averag e 
molecular gas density is ~10® cm“^ ( Scoville et al.|2014 ), 
comparable to the conditions in a dense star-tormiiig core 
inside a giant molecular cloud. Indeed, Ar p 220 is un¬ 
dergoi ng very high rates of star formation. [Smith et al. 
(1998) estimate the circumnuclear star formation rate at 
50-800 Mf7> yr~^ u sing the radio supernovae rate, while 


Farrah et al.| (2003) use the far-infrared luminosity to es¬ 
timate a global star formation rate of 240±30 Mq yr“^. 

Outside of the heavily extincted region near the nuclei, 
star formation activ ity was observed with HST ACS. 
Wilson et al. (2006) found two populations of massive 
star clusters: young (<10 Myr) massive clusters within 
2.3 kpc of the nuclei and intermediate age (~I00 Myr) 
clusters beyond this radius. 

In this paper, we investigate the connection between 
the nuclear energy sources and the large-scale outflow. 
We use HST/WFC3 with narrowband filters to create the 
first high spatial resolution optical emission line maps of 
Arp 220. We present maps of [OIII], H/3, and Ha-|-[NH] 
to resolve gas motions in Arp 220. We discover a bubble 
feature in Hq!-|-[NH] (section]^ and discuss its origins in 
section ID 

Throughout this work, we assume a distance of 77 Mpc, 
which gives 1”=375 pc. 

2. OBSERVATIONS 
2.1. HST 
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Table 1 

HST WFC3-UVIS Exposures 


Filter 

Bandpass 

Central A (A) Width (A) 

Line 

Exposures 

Continuum 
Eilter Scale*^ 

Flux Error 

1(T^ (erg s~^ cm~^ A“^ arcsec“^) 

FQ492N 

4933 

114 

H/3 

2x1090 s, 6x920 s 

F467M 

0.42 

3.9 X 10-1® 

FQ508N 

5091 

131 

[OIII] 4959, 5007 

2x 1050 s, 6x920 s 

F467M 

0.46 

3.8 X 10-1® 

F665N 

6656 

131 

Ha -1- [Nil] 6548, 6583 

2x200 s 

F621M 

0.16 

1.4 X 10-11’ 

F680N 

6877 

371 

[SII]'^ 6717, 6731 

2x525 s 

F621M 

0.66 

7.7 X 10-1® 

F467M 

4683 

201 

blue continuum 

2x775 s 




F621M 

6219 

609 

red continuum 

2x155 s 





^ Scaling factor of the continuum used for continuum subtraction. 

^ Flux errors are estimated based on the dispersion of the sky background in each image, as measured on a section of sky away from the galaxy. 
^ Not used for analysis due to shallowness of image; see Appendix 


The optical imaging data of Arp 220 used in this work 
are new observations obtained with HST Wide Field 
Camera 3 (WFC3-UVIS) on 2012 September 22 (GO- 
12552, PI: Lisa Kewley). We selected filters whose band- 
passes each cover one or more of several strong optical 
emission lines commonly used for calculating gas-phase 
metallicity and shock excitation: H/3 (filter FQ492N); 
[OIIIJAA 4959, 5007 (rest) (FQ508N); Ha -b [NH]AA 6548, 
6583 (F665N); and [SH]AA 6717, 6731 (F680N). Blue 
(F467M) and red (F621M) continuum exposures were 
also taken to allow calculation of the line flux across the 
extent of the system. Figure 1^ shows the filter band- 
passes superimposed over the SDSS 7 spectrum of Arp 
220. The filter and exposure data are summarized in Ta¬ 
ble [l] along with the continuum bandpass used with each 
emission line imag e and the la R MS errors (discussed 
iEO’ 

3t viei 


2.1.2) With a pixel size of 


further in sections |2.1.1 and 

0.0396”, the field of view (FOV) of WFC3-UVIS is 164” 
X 176” as measured on the image. The FQ, or quad, 
filters, each occupy one-quarter of the FOV, so the FOV 
for these filters is approximately 82” x 88”. 



4500 5000 5500 6000 6500 7000 7500 
Wavelength (A) 


Figure 2. SDSS 7 spectrum of Arp 220 with the WFC3-UVIS 
filter bandpasses overlaid. The left axis shows the throughput of 
each WFC3 filter (colors are as shown in the legend) while the right 
axis shows the flux of the SDSS spectrum (in black). 


Data reduction was performed using the STScI 
pipeline, which utilizes “On-The-Fly Reprocessing” 
(OTFR). This processes the raw data files using the best 
calibration files and software available at the time of re¬ 
trieval from the archive and performs basic bias and dark 


subtraction and flat fielding. It also flags bad pixels and 
performs initial cosmic ray removal. The data were all 
retrieved from the archive in late 2012. 

The reduce d dithered exposu res were combined using 
AstroDrizzle (Avila et al. 2012). AstroDrizzle, the suc¬ 
cessor to MultiDrizzle, provides distortion corrections 
and combines multiple images using a technique which 
allows for image subsampling. Sky subtraction was per¬ 
formed using the mode of the image, in order to not bias 
the sky value by the extended nebulosity of the galaxy. 
Charge transfer efficiency (GTE) trails caused by ineffi¬ 
cient electron shuffling during readout were evident on 
the initial combined images. Due to the high counts in¬ 
volved, they were especially problematic in the pixels 
along which cosmic ray hits were shuffled. The charge in 
the first smeared pixel amounted to roughly 1% of the 
initial high pixel value and dropped off rapidly as the 
charge was further shuffled along the detector. However, 
the GTE trail left by each cosmic ray was up to 20 pix¬ 
els long for the strongest cosmic rays. These pixels were 
masked in AstroDrizzle, as the data were downloaded 
and reduced before the GTE correction software from 
STScI became available. All other AstroDrizzle param¬ 
eters were set to the default values. The final combined 
images were registered to a common reference frame us¬ 
ing the AstroDrizzle tasks tweakreg and tweakback and 
redrizzled onto the same pixel grid, to allow for image 
combination. 

The EQ492N and FQ508N filters, covering H/3 and 
[OHI], had a sufficient number of exposures to allow 
for cosmic ray removal (see table . The other filters 
(F467M, F621M, F665N, F680N) nad two dithered ex¬ 
posures each; the cosmic ray cleaning provided by As¬ 
troDrizzle for these was not robust, especially in the gap 
between the two detector chips. The final combined im¬ 
ages for these filters were further cleaned of cosmic rays 
using a Laplacian edge detection algorithm, which iden- 
tifies cosmic rays ba sed on the sharpness of their edges 
(van Dokkum 2001). This method, which removes only 
the positive noise spikes, can alter the noise distribution 
of the image. The resulting noise distribution may be 
asymmetric and skewed towards negative values, which 
can be problematic for low signal-to-noise ratio (S/N) re¬ 
gions. We examined the cleaned images and found that 
only the F665N image is affected. Because we focus on 
the high S/N regions in this filter, this does not affect 
our results. 

The bandpass of one of the continuum images, through 
filter F621M, contains [OI]A 6300A, a strong emission 
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line in shock-excited gas. The strength of this line varies 
with metallicity and excitation source and thus could 
contaminate calculations of emission line ratios. The 
contribution of the flux of this line to the total flux in 
the continuum image was estimated using the Sloan Dig¬ 
ital Sky Survey (SDSS) spectrum of Arp 220. The SDSS 
fiber covers the inner 3” of the galaxy, or roughly 1 kpc. 
The ratio of flux in the emission line to flux in the contin¬ 
uum contained in the bandpass was roughly 1%, or less 
than the average error in the images. Thus, the contri¬ 
bution of the line to the continuum image was taken to 
be negligible and no correction to the continuum image 
was made. 

The data were binned in square bins 3 pixels wide to 
increase the S/N. At th e redshift of 0.018126±0.000023 
(de Vaucouleurs et al. 1991), the WFC3 pixel scale of 
0.039(5 "/pixel corresponds to 15 pc. Extragalactic giant 


HII regions range from 50 to 300 pc in diameter (Oey 


et al. 2003). The binned resolution elements measure 
45 pc across and thus allow full sampling of all but the 
smallest HII regions while still allowing recovery of low 
signal. 


2.1.1. Continuum Subtraction 


One of the most sensitive calibration tasks is subtrac¬ 
tion of the stellar continuum from the emission line im¬ 
ages. Small changes in the factor used to scale the con¬ 
tinuum before subtracting it from the emission line image 
propagate into a change in the flux of the final emission 
line image. This can have a large effect on the emission 
line flux and ratio maps. 

We determined the optimal scale factor between the 
continuum image and the emis sion line image us ing the 
“skewness transition method” (Hong et al. 2014). This 
method uses the symmetry of the flux distributioh within 
an aperture to identify the ideal scaling factor that nei¬ 
ther under nor over subtracts the continuum. It accounts 
for spatially dependent luminosity changes, bnt not spa¬ 
tially dependent color changes. The location of the aper¬ 
ture used to calculate the scaling factor influences the 
derived value, and allows us to estimate an error on the 
scaling factor. The error on the scaling factor can be 
taken as 0.05 for each filter. 

We verified the accuracy of this method via physical 
arguments. For a star cluster with zero extinction, and 
assuming Case B recombination, the ratio of Ha to H/? 
is 2.86. However, Arp 220 has high levels of patchy ex¬ 
tinction and a high and varying ratio of [Nil] to Ha, so 
we would not expect this ratio to hold across most of 
this system. We are able to check the ratio between the 
two images using a star cluster with significant Ha-|- [NH] 
and H/3 emissi on, and a measured E(B-V) of 0 (Wil¬ 


son et al. 
ratio of 


2006, their cluster 62). We assumed that the 
SI III to Ha is roug hly that for an HII region 


(^0.3, [Kewley et al. 2006a) and adjusted the scaling 
factor for the continuum subtraction from the H/3 im¬ 
age until we achieve a ratio equal to the Balmer decre¬ 
ment. With this method, we calculate a new scaling fac¬ 
tor of 0.49±0.02. This scaling factor is generally consis¬ 
tent with that obtained through the skewness transition 
method (0.42±0.05), and so we adopt the scaling factors 
from the skewness transition method. Our scaling factors 
are given in table 


2.1.2. Error estimation 

Flux errors are estimated based on the standard devi¬ 
ation of the background in each image. The dispersion 
varies slightly across the field, so to obtain a representa¬ 
tive value this was measured in a region centered within 
2” of the detected emission line flux in each filter, but 
not overlapping it. We show the Icr errors in Table 
We adopt a S/N cutoff of Aa in all emission line images. 
However, the F680N image was too shallow and we were 
unable to recover [SH] flux above the 4cr limit (see Ap¬ 
pendix). We thus drop this image from further analysis. 


2.2. Other data 
2.2.1. H 2 data 


Near-infrared (H and K band) integral field spec¬ 
troscopy (IFS) data were taken 2007 March 7 and April 
18-21 with SINFONI, a NIR integral field spectrometer 


on th e Very Large Telescope (VLT) (Eisenhauer et al. 


20031. The curvature-based adaptive optics system was 


operated in laser guide star mode, though without tip-tilt 
correction as no suitable tip-tilt star was available. This 
gave a nearly symmetric Gaussian point spread function 
of 0.30” X 0.31”. The 0.05” xO.10” plate scale was used, 
giving a FOV of 3.2” x 3.2”. Further details on the ob- 
servations, data re duction, and line fitting are given in 


Engel et al. (2011). 


2.2.2. NICMOS 

Near-infrared imaging using HST/NICMOS was ob¬ 
tained from the HST archives. The F160W filter on cam¬ 
era NIC2 was used to image Arp 220 in the H-band on 
10 January 2004 (proposal 97 26, PI: R. Ma iolino) for a 
total exposure time of 600 s (Cresci et al.|2007|. The 
F222M filter (K band) on NI(J2 was used to image the 
sy stem on 4 April 1997 for a total exposure time of 1024 
s (Scoville et al. ||1998l proposal 7116, PI: N. Scoville). 
The pixel scale on the NIC2 camera is 0.076” while the 
FOV is 19.2”. 


2.3. Registration 


There are no stars suitable for image registration in the 
field of view (FOV ) . Inst ead, bright super star clusters 
from IWilson et al.| (|2006|) were used to align the WFC3 
and iNlOMOS datasets. Three bright, well-separated 
clusters were used, allowing a correction for both shift 
and rotation between the images. The absolute astromet¬ 
ric registration was set by the coo r dinate s of the west¬ 
ern nuclei given by Scoville et al. (1998), precessed to 
J2000.0. 


The SINFONI data have a smaller FOV, limited to 
the circumnuclear region. The optical super star clusters 
(SSCs) previously used for alignment were out of this 
FOV. Instead, the IFU data were aligned directly with 
the NICMOS F222M image using the nuclei and another 
bright infrared source. 


3. MORPHOLOGY 


3.1. Optical emission line structure 

We present emission line maps of [OHI], H/3, and 
Ha-|-[NH] in figure]^ Each image is 45” x 45”. A 4cr S/N 
cut has been appliS, and the logarithm of the emission 
is displayed in order to show the full range of emission. 
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Figure 3. Maps of the log([OIII]) (left), log{H/3) (middle), and log(Ha+[NII]) (right) emission across the galaxy, rotated to acc ount for 
a 77° P.A. s o that north is up and east is left. The locations of the eastern and western nuclei, as determined in the near IR by |Scoville| 
|et al.| ( |1998| l and processed to J2000.0, are shown as red diamonds. A box, showing the size of the zoomed figure]^ is shown on the right 
panel. A 4cr S/N cut has been applied. 
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Figure 4. Map of the log([OIII]) (left), log(H/3) (middle), and log(HQ:-|-[NII]) (right) emission near the nuclear region, zoomed in on a 
central 6” x 6” region. As in the full map, a 4cr S/N cut has been applied, the locations of the eastern and western nuclei are shown as 
red diamonds, and the maps have been rotated so that north is up and east is left. 


The maps have been rotated to account for a 77° P.A., 
so that north is up a nd east is left. The NIR locations 
of the double nuclei (|Scoville et al.||1998|) are shown as 
red diamonds. Some residual cosmic rays are visible in 
the chip gap of the Hq!-|-[NII] (right) panel figure]^ as 
detailed in section [21 

The imaging reveals widespread [OIII] and emis¬ 
sion across the central disk, roughly 8 kpc in diameter. 
Faint tidal tails are visible to the northeast. A dust lane 
runs from the northeast to the southwest, obscuring the 
double nuclei of the system, which are only visible in the 
X-ray, IR, and radio. Though visible across the main 
disk of the system, the [OIII] and H/3 emission is con¬ 
centrated near the optical center, which is roughly 1.5” 
offset to the northwest from the western nucleus. The 
Ha-|-[NII] emission is much more tightly concentrated 
near the optical center. We also present emission line 
maps zoomed into the nuclear region in figured 

The most striking feature in the emission Ime maps 
is a ring- or bubble-shaped structure, seen in Hq;-|-[NII] 
emission (figure right panel). It is unclear whether 
the feature is due to strong emission in [Nil], Ha, or 
both. Strong Ha emission indicates high rates of star 
formation, while high [Nil] emis sion generally indicat es 
the presence of shocked gas (see[Kewley et al.[ ([2006a[)). 
The ring measures 1.6” in diairieter, or 601) pc at the 
distance of Arp 220, with a P.A. of 316° as measured 
from the western nucleus to the bubble center. There 
appears to be a faint smaller ring, inside the larger one, 
on the side closest to the nuclei. This is the first time the 


Ha and [Nil] peaks reported in Arribas et al. (2001 their 
figure 1, approximately corresponding to their hbers 110- 
114) have been resolved into a bubble morphology. This 
bubble fe ature and its environment are discussed further 
in section [3]3l 


3.2. Emission line ratios 


The ratios of strong optical emission lines are com¬ 
monly used as a diagnosti c to determine the excitation 
mechanism of ionized gas ([Baldwin et al. 1981[ [Kewley 


et al.|200T Kauffmann et al.|2003 ' Kewley etm. [20066 [T 
photoionization by massive stars, excitation via shocks, 
or ionization from an AGN. The most widely used diag¬ 
nostics require the ratio of both [OIII] to H/? and either 


[Nil], [SH], or [01]A 6300A to Ha. 

In figure we present maps of log([OIII]/H/3), which 
were createa from the left and center panels of figure [^ 
A S/N cut of 4cr was first applied to each individual emis¬ 
sion line map and the resulting ratio map was smoothed 
using a median filter (kernel=3) in order to empha¬ 
size larger-scale features. The right panel of figure [^ is 
the same, but cropped to more closely show the region 
around the bubble feature. Ha-|-[NH] contours are over¬ 
laid in black on the right panel, to show the location of 
the bubble feature. 

Some structure is visible in the smoothed maps, and 
is more readily apparent in the normalized histograms of 
log( [OIII]/H/3) shown in hgure 1^ Values for the entire 
system are shown in blue, while values within just the 
Ha-|-[NH] ring feature are shown in red. Binned regions 
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Figure 5. Median smoothed ma ps (kernel=3) of log([OIII]/H/3). Red diamonds show the positions of the nuclei, identified in the NIR 
and given in|Scoville et al.| (|1998|l. Left: Full view of system. Right: Same as left, but zoomed into the optical center of the system. In 
addition, contours showing the bubble feature, in Ha+[NII] emission, are shown in black. North is up, east is left in both panels. 


are considered part of the ring feature if they fall within 
1” of the center of the ring. The mean log([OIII]/H,d) 
for the entire system is 0.05±0.11, or 0.22±0.07 for the 
bubble region (shown in blue and red, respectively, above 
the histograms). The log([OIII]/H/3) values in the ring 
feature appear to be slightly higher than in the system 
as a whole, though the difference is not significant at the 
Icr level. 

The value of log([OIII]/H/3) will be used in conjunction 
with literature values of log([NII ] /Ha) to analyze the gas 
excitation mechanism in section T4.2.1l 



log([0lll]/H/S) 


Figure 6. Normalized histograms of log([OIII]/H/3 values of the 
map of Arp 220, so that the area under each histogram is 1. The 
histogram in blue shows all log([OIII]/H/3) values across the sys¬ 
tem (N=18676), while the histogram in red shows only the values 
within the Ha-|-[NII] bubble feature region (N=236). The mean 
and standard deviation of each histogram are shown at top. 


3.3. The bubble in a multi-wavelength context 

In this section, we present the first high spatial res¬ 
olution multi-wavelength maps of Arp 220 combining 
optical, NIR, and X-ray images. A 3-color image, 
shown in figure was made using the optical continuum 
(HST/WFC3 F621M) image, Ha-|-[NII] emission image, 


and the H-band continuum (NICMOS FI60W) image. In 
this image, the blue shows the optical stellar continuum, 
subject to extinction from dust in the system, especially 
in the circumnuclear region; the green shows the sum of 
the Ha and [Nil] emission lines, which can be excited 
by recent star formation, shocked gas, or AGN excita¬ 
tion; and the red shows a combination of the near-IR 
continuum, which is less affected by dust extinction. 

The X-ray structure seen on ~10 kpc scales, which 
traces the hot gas outflow, appears to be roughly aligned 
with the P.A. of the bubbl e feature with respect to the 


western nucleus (figure 
2) re ports a P.A. of 3 


while 

331 


McDowell et al. 


T \ . [Heckman et al.| (|i996| figure 
5° fo r the northwestern plume, 


2003 figure 1) report a P.A. of 


or the same feature. However, the structure is 


wide and faint enough that determining an ex act P.A. 
is difficult. In addition, Heckman et al. (1996) report a 
large-scale X-ray structure with a P.A. of 285°, which 
is inconsistent with the P.A. of the bubble feature. In 
figure we show the mid-scale X-ray structure as re- 
porteclDy Heckman et al. (1996). The binned pixel size 
in the X-ray data (1-5") is too large to resolve whether 
the outflow is coming from one of the two nuclei or the 
kpc nuclear gas disk. Higher spatial resolution obser¬ 
vations reveal a hard X-ray source roughly aligned with 
the western nucleus (Clements et al. 2002, source X-1). 
Therefore, we proceed with the assumption that the X- 
ray outflows originate in the western nucleus. 

The southeastern side of the bubble appears to curve 
around the western nucleus. This curve is echoed by the 
NIR emission surrounding the wester n nucleus, which 


likely traces the circumnuclear disk (Sakamoto et al. 


1999). While little optical emission is seen near the nu¬ 


clei, due to high dust extinction, there are peaks in the 
optical emission lines (Ha-I-[Nil] and [OHI]) roughly co¬ 
incident with the NIR nuclei. As the images were regis¬ 
tered based on the coordinates of SSCs throughout the 
galaxy, not the nuclei, this alignment is remarkable. 

A gap in the ring feature is evident on the side away 
from the nuclei, which could be due to either high extinc¬ 
tion or an actual gap in emission. If the gap is caused 
by extinction, and the obscuring material contains hot 
dust, we would expect to see more NIR emission coinci- 
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Figure 7. 3 color image of Arp 220. The blue image shows the optical continuum (HST/WFC3 F621M), the green shows Ha+[NII] 

(continuum-subtracted HST/WFC3 F665N), and the red shows the H-band infrared continuum (HST/NICMOS F160W). The major axis 
of the soft X-ray structure (white vector) and locations of the super star clusters (SSCs, white crosses) are indicated. The left panel shows 
the main disk of the system. The scaling of the red image is adjusted to emphasize the faint emission; the chip gap is visible down the left 
side of the image. The right image shows the nuclear region, including the double nuclei and the Hq-|-[NII] bubble feature. The peaks in 
the infrared mark the positions of the two nuclei, and are shown as white circles. The scaling of the red image has been adjusted to show 
the structure of the strong emission near the nuclei. 


dent with the gap. The NIR emission does not peak in 
this area, indicating that this hole in emission is not due 
to extinction from material containing hot dust. In addi¬ 
tion, the mid-scale northwestern X-ray plume is aligned 
with this gap, suggesting that a stream of hot X-ray emit¬ 
ting gas is clearing a hole in the optical emission line gas, 
disrupting the process that formed these lines. This hot 
gas extends for 10-15 kpc on each side of the nucleus 
(McDowell et al.||2003). 

INo SSCs appear to be coincident with the bubble fea¬ 
ture. Figure shows the locations of the optical and in¬ 
frared selects SSC s as white crosses ( Wilson et al.|200j^ 


Scoville et al.|1998). We discuss this further m section |4 


3 . 4 . Kinematics from the NIR 


Engel et al. (2011) studied the NIR kinematics of the 
circumnuclear region (r~1.5”, centered on the western 
nucleus) using the H 2 1-0 S(l) line. The velocity struc¬ 
ture presented there shows a clear rotational signature, 
with the rotational axis roughly in the NW/SE direction. 
The NIR velocity dispersion map is shown in figure 
with the optical Ha-(-[NII] flux contours ov erlaid. This 
map d iffers slightly from that presented by [Engel et al.| 
(2011) in that the instrumental broadening has been cor¬ 


rected for here; thus, the dispersion values are slightly 
lower than those presented earlier. The velocity disper¬ 
sion of the H 2 is enhanced near the bubble, and reaches 
a maximum of 250 km s“^ at the base of the western 
arm. The velocity map shows no evidence for a systemic 
velocity shift at this location, though the rotational axis 
is not perfectly aligned with the region of enhanced ve¬ 
locity dispersion. This may indicate that the higher ve- 



Figure 8. NIR velocity dispersion, calculated from the H 2 1- 
0 S(l) line and corrected for instrumental broadening. Contours 
show the 7 (t Ha+[NII] flux, smoothed for display purposes. The 
nuclei are marked as cyan diamonds. North is up and east is left. 


locity dispersion is linked with the bubble feature. If we 
assume an outflow is launched from a source at the cen¬ 
ter of the bubble and that the outflow velocity is equal 
to the peak of the dispersion, then the age of the bubble 
is ~2.3 Myr. If the powering source is in the nucleus, a 
kinematic age cannot easily be determined. The current 
kinematic maps do not cover the bubble’s full spatial ex¬ 
tent. Additional data are necessary to characterize the 
dynamics of the region. 

3 . 5 . Bubble energetics 





















We estimate the mass and energy in the bubble using 
and assuming Case B recombination and T=10'* 
K (Osterbrock 19891. Using the continuum-subtracted 
HQ!-|-[iNllJ image, we find a total integrated flux in the 
bubble region of 2.6x10“^"^ erg s“^ cm“^. With the Ha 


to Ha-|-[NII] ratio from the bubble region of the Arribas 


et al. 


(20011 IFS study (/Jq^P[%77j= 0.37), we calculate an 

ionized gas mass of 2.5x10^ Mq 
The Ha flux used here has not yet been corrected for 
extinction. We can derive the extinction in the bubble 
region via the Balmer decrement, by taking the ratio of 
the Ha to H,d maps and comparing this value to 2.86, 
the theoretical ratio under Case B recombination. This 
calculation is very sen sitive to the continuum subtrac¬ 
tion for H/3 (see section 2.1.1), but we can derive a range 
of acceptable extinction values based on the error in the 
continuum subtraction scaling factor. The mean con¬ 
tinuum scaling factor for H/3, 0.42, gives E(B-V)=0.11, 
though a range in E(B-V) of 0.00-0.48 is also consistent. 
Given the mean extinction, we derive a dereddened ion¬ 
ized gas mass of Mio„=3.2x 10^ Mq nj^, with a range of 
2.5-7.8x10^ Mq We adopt the mass derived from 
the mean value of the H/3 continuum scaling factor in the 
following discussion. 

The [SH] doublet ratio 6317A/633lA can be used to 
constrain the electron density of the ionized gas. The 
MPA-JHU SDSS DR7 data product catalog gives a ra- 
tio of 1.46, which places Ug in the lower density limit 
(Osterbrock 1989), or less than ^10^ c m~^. Previous 
work on the ionized bubble in NGC 3079 (Veilleux et al. 


1994 Gecil et al. 2001) constrained ne in that system 


Rupke & Veilleux (2013b I adopt ne=10 


to 5-100 cm^ 

cm“^ for their analysis of the ionized winds in a sample of 
ULlRGs, consistent with the NGG 3079 limits. We adopt 
ne=10 cm“^ here to facilitate comparison with the lit¬ 
erature. Assuming t he upper limit of n p =10 cm ~^ gives 


M,o„=3.2xl0® Mq. iRupke & Veilleuxl (|2013bl Table 6) 


find ionized gas masses in their sample of ULIRG winds 
in the range iqS-^^-s.o? so the Arp 220 bubble mass 
is similar to these other ULIRG outflows. 

We estimate the kinetic energy in an expanding shell 
by assuming that the observed velocity dispersion in H 2 
represents the shell velocity, or Vsheii=‘^50 km s“^. We 
find a kinetic energy of 2.Ox 10^^ erg. Given a dynamical 
age of 2.3 Myr (section 3.4) and assuming constant en¬ 


ergy injection, we calculate an energy injection rate from 
the ionized gas mass of 2.8x10^° erg s“^. 

We are only observing the ionized gas and the neutral 
gas may be a substantial fraction of the total gas mass 
in the bubble feature. To include a correction for the 


neutral hyd rogen, we used the measurements of Mundell 
et al. (2001), who measured hydrogen absorption against 
the continuum in 1.4 GHz observations. The structure 
they identify as spur T is coincident with the bubble fea¬ 
ture, though they are unable to measure significant H 
absorption in this region. They place an upper limit of 
Nh < 1.3x10^° Ts(K) cm“^. With their assumption 
that Tg = 100 K, that gives an upper limit of Njj < 
1.3 xlO^^ cm“^, and an upper limit on the total neu¬ 
tral hydrogen mass in the bubble region of 1.2x10® Mq. 
Using the same velocity assumptions as for the ionized 


® http://www.mpa-garching.mpg.de/SDSS/DR7/ 



Figur e 9. Left: Last panel of figure 2 from |Sutherland fc Bicknell| 
l|2007|, reproduced here with permission. These simulations model 
the iriteraction of a supersonic jet with an inhomogenous ISM. 
The figure shows the 95 kyr time-step of the mid-plane density 
slices of the simulation, with darker colors representing less dense 
material. The panel is 1 kpc on a side. Right: 3-color image of 
the bubble feature, as in the right panel of figure rotated and 
zoomed for comparison with the simulation. The western nucleus 
is centered on the left side of the panel and the opening in the 
bubble is centered on the right side of the panel. This panel is also 
1 kpc on a side. 


gas, the upper limit of kinetic energy in the neutral gas is 
3.8x 10®® erg, and the upper limit on the energy injection 
rate is 5.3x10^^ erg s“^. 

The total neutral-fionized mass in the bubble fea¬ 
ture is thus M<1.2xl0® Mq, the total kinetic energy 
is KE<4.0xl0®® erg, and the energy injection rate is 
dE/dt~5.6xlO^^ erg s“^. We compare these rates to 
possible sources in the following section. 

4. ORIGIN OF THE BUBBLE 

Given the bubble’s position and morphology, there are 
two possible points of origin. First, it could be the prod¬ 
uct of a jet or outflow originating in one of the two nuclei 
or the molecular disk immediately surrounding them. Al¬ 
ternatively, the bubble could originate from a source at 
its center. Either a massive star cluster or a microquasar 
(i.e. an ultraluminous X-ray source, ULX) could drive a 
traditional expanding superbubble. Given that the nu¬ 
clei are deeply embedded in dense gas and are more dis¬ 
tant, the energies necessary to drive an outflow from the 
nuclear disk versus from the bubble’s center are v ery dif- 
ferent. We lay out the possibilities in sections [4T]and|4.2[ 
and discuss and compare them further in section |4.3[ 

4.1. Nuclear origins 

The first possibility is that the bubble is being driven 
from a nucleus or the circumnuclear region, either from 
an AGN(s) or from a strong nuclear starburst. 


4.1.1. Jet or outflow from an AGN(s) 

The observed morphology of the bubble feature 
matches that produced in s ome AGN jet simulations. 
Sutherland & Bicknell (2007) modeled the interaction of 
a light supersonic jet with the inhomogeneous and tur- 
bulently supported ISM in an early-type galaxy. Their 
simulated ISM followed a lognormal single-point density 
distribution and a power-law density structure, with an 
average jet/ISM number density ratio of ~10“^. The jet 
had a kinetic luminosity of ~10^® erg s“^, and radiative 
cooling was implemente d via the MAPPINGS phot oion- 
ization and shock code (Sutherland fc Dopita|[l9M ). 
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The evolution of the interaction goes through several 
stages, including the development of a shocked bubble 
as the jet interacts with the ISM followed by a break¬ 
through of the jet through the outer boundary of the 
bubble. The simulation mid-plane density slices post¬ 
breakout (their figure 2, and reproduced here, in figure]^ 
resemble the Arp 220 bubble. However, because this is 
a generic jet/ISM simulation, this is an idealized com¬ 
parison; in particular, the ISM of Arp 220 is likely to be 
denser than that assumed here. Nonetheless, the mor¬ 
phological similarities suggest the possible presence of 
an AGN and its associated outflow. 

The proximity of the nuclei to each other (360 pc) 
suggests that they should be subject to strong gravita- 


AGN (Rupke et al. 2005a Tremonti et al. 2007), but 
may b e consistent with velocitie s from lowduminosity 


an AGN has been sparse. Previous X-ray 

studies (Rieke 

1988 

Dermer et al. 1119971 llwasawall 19991 llwasawa et al. 

2001 

Clements et al. 2002) ruled out a 

bright hard X- 


must be highly obscured, with a lower limit on the nu¬ 
clear column density of 10^® cm“^. However, evidence for 
at least one highly obscured AGN in Arp 220 has been 


mounting in recent years. A maser (Aalto et al. 


and a rotating massive molecular disk (UownesTz 


2007) have been found, and Contini (2013) find evidence 


2009 


Ickar 


of an AGN in th eir modeling of the nuclear optical and 
infrared spectra. Paggi et al. (2013) find strong Fe-Ka 
line emission in the vicinity of the two nuclei, concen¬ 
trated most strongly on the western nucleus. From their 
broadband X-ray data, they estimate lower limits on the 
AGN bolometric luminosities of 5.2x10^^ erg s“^ (west¬ 
ern nucleus) and 2x1 0^^ erg s~^ (eastern nucleus). New 


ALMA observations (Wilson et al. 2014) show a lumi¬ 
nosity surface density in excess of KT^L© kpc“^ at the 
western nucleus, sufficiently bright to require the pres¬ 
ence of a hot starburst, a Co mpton thick AGN, or both. 

As shown in section |3.5[ the kinetic energy input 
needed to form the bubble is dE/dt=5.6x 10^^ erg s“^ 
(upper limit, ionized-)-neutral medium). Taking this as 
a ratio with respect to the bolometric luminosity of the 
suspected AGN in the western nucleus gives an energy 
deposition rate of dE/dt / 1->agn = 1.1%, whic h is an or¬ 


der of magnitu de higher than other observed (Rupke & 
WiUe^ 2013b|) or mod eled AGN energy feedback rates 


j Hopkins fe Elvis||2010| . However, these models are de- 
pendent on the traction of gas in the hot ISM and are 
uncertain enough to encompass our observed value. 

The evidence for an obscured AGN driving the bubble 
is favorable but not absolute. The morphological simi¬ 
larities to models and the energetics lend support to the 
theory that this bubble is formed by the interaction of 
a jet or outflow from an AGN with the galaxy’s ISM. 
This is further supported by the large-scale X-ray struc¬ 
ture of the galaxy, which is remarkably well aligned with 
not just the bubble feature, but also the gap in the side 
of the bubble away from the nuclei (figure Q . However, 
the western nucleus (the most likely AGN host) has a gas 
accretion disk that is misaligned with the supposed direc¬ 
tion of the jet and bubble. This is not definitive evidence 
against an AGN outflow, because the observed gas disk 
has a scale of 100 pc and any inner accretion disk may 
have a different orientation. Additionally, the outflow 
velocities in this system are on the low end for typical 


AGN (Christopo ulou et al. 1997 Kewley et al. 2006b 
Hicks et al.||2U13). We view an AGiN jet or outflow as a 


plausible hypothesis for this bubble’s origin. 

4.1.2. Collimated outflow from a nuclear starburst 

The prevailing explanation for the large-scale super¬ 
wind structure as seen in X-ra ys is a massive starburs t 
concentrated near the nuclei (Heckman et al. 1996). 


Smith et al. (1998) used VLB! continuum imaging to 
resolve a strong, compact OH maser near the western 
nucleus into a number of radio supernovae within a 75 
X 150 pc region. Using the lifetime of these objects, 
they estimate a nuclear star formation rate of 50-800 Mq 
yr“^. Such high levels of star formation in the confines 
of the nuclear gas disk could potentially drive a starburst 
wind. If the starburst wind is able to break through the 
gas disk, the dense gas in the disk will further direct the 
outflow along the maximum pressure gradient. 

In section [STsl we showed that the kinetic energy input 
needed to form the bubble is dE/dt=5.6x 10^^ erg s“^ 
(upper limit, ionized-)-neutral medium). If we assume 
a SFR of 240 Mq yr“^ (equal to the global SFR and 
well within the nuclear SFR range), we can compare the 
energy in the wind to the m echanical luminosity in th e 
nuclear starburst. Following Rupke & Veilleux (2013b), 
we calculate the mechanical luminosity in arontmuous 
starbu rst using Starburst99 models from [Leitherer et al.| 
(1999). We apply a correction for a Salpeter IMF, with 
a lower mass cut at O.I M©. Thus, the energy deposition 
rate is dE/dt / dE/dtssgg = 0.8% (upper limit). This is 
in line with that observed in other starburst-dominated 


ULIRGs (3%-I3%, Rupke & Veilleux 20I3b). It is thus 
also possible that the bubble feature is driven by the 
strong nuclear starburst. 


4.2. Off-nuclear origins 

The bubble may originate from a source at its center 
rather than from the nucleus. It could be the product 
of a massive star cluster, or of a microquasar or ULX 
object. 

4.2.1. Massive star cluster 

Structures morphologically similar to the Arp 220 bub¬ 
ble feature, known as superbubbles, are believed to be 
formed by supernovae winds from a m assive star cluste r 
pushing out into the surrounding ISM (Cash et al.|1980 ). 
The Arp 220 bubble feature demonstrates some evidence 
consistent with a classical superbubble, but lacks evi¬ 
dence of a sufficiently massive star cluster driving the 
outflow. 

The Arp 220 bubble shows evidence of shock excita¬ 
tion, which is consistent with a supernovae-driven out¬ 
flow from a massive young cluster. Such outflows shock 
as they encounter the surrounding ISM and produce an 
increasing ratio of [Nil] to Ha emission at a fixed value of 
log([O HI]/H/I) (|Baldwin et al.||I981[ |Kewley et al.||2001 
2006bl) . The value of log(|Uili|/H/:/) is fairly constant 
acro ss t he center of the galaxy at 0.22±0.07 (see sec¬ 
tion 3.2 0.05±0.11 for the galaxy as a whole). The value 
of log([NH]/Ha) at the maximum starburst line, or the 
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value at which shock heat ing is thought to begin to dom¬ 
inate (Kewley et al.|2001 1, is -0.16±0.04 for this value of 
log([ Olll|/H;d) (though there are metallicity effects; see, 
e.g., IDopita et al.1|2000[ IDenicolo et al.||2002| [Kewley ^ 
Dopita||2002| [Kewley et ai.||2006bp. Though we cannot 
separate the [1N 11| and tia em i ssion to calculate a ratio 
with our data, Arribas et al. (2001) show optical IFU 
spectroscopy data taken with the INTEGRAL instru¬ 
ment on the William Herschel telescope. Their data are 
lower spatial resolution (fiber diameter^O.9”, R~2000) 
and they identify three gas components via spectral line 
profile decomposition. Spatially, our bubble feature is as¬ 
sociated with their component B. They find an average 
value of log([Nll]/Ha) of 0.48±0.17 for this component. 
This is consistent with shocked gas in a shell surrounding 
a massive star cluster that is driving the outflow. 

Although our feature resembles typical massive star 
cluster superbubbles, no massive star clusters are visible 
in the center of the bubble (figure [7). Examination of 
both the optical continuum image (\WC3-UVIS F621M) 
and the H-band image (NICMOS F160W) reveals that 
the flux interior to the ring feature is high enough that a 
10® Mq cluster would be hidden by the background. To 
determine this limit, we measured the flux of other 10® 
Mf7^ c l usters outside the bubble as identihed by [Wilson 


et al. (2006 , clusters 2, 3, 4, 5, 11, and 12). The mean 


-2 


flux for a 10® Mq cluster is 3.52±1.50 x 10 ergs cm 
s-i A-i for WFC3 F621M and 1.77±0.86 xlO"!® ergs 
cm“^ s“^ for NICMOS F160W, using an aperture 
of radius 0.17”. 

Analysis of the energetics provides additional evidence 
against an outflow driven by a massive young cluster 
at the bubble center. The mechanical l uminosity pro¬ 


duced by a 10® Mq cluster (^10"^® erg s ^ Leitherer et al. 


1999) is a factor of ~100 below the energy injection rate 
(~10^^ erg s“^) we derived for the bubble. Thus, while 
we cannot rule out the existence of an undiscovered clus¬ 
ter smaller than 10® Mq , it is unlikely that it is the origin 
of the bubble feature. 


4 . 2 . 2 . ULX bubble 

Ultraluminous X-ray sources (ULX) are compact X- 
ray sources whose luminosities exceed the Eddington 
luminosity of an isotropically-emitting stellar-mass ob¬ 
ject. By definition, they are located away from the nu¬ 
cleus of their host galaxy. The central excitation source 
is still unknown, though their observed variability has 
been used to hypothesize these as accreting stellar- or 
intermediate-mass black holes (Colbert & Mushotzky 
1999 Kubota et al. 2001). They can exhibit either soft 
or hard X-ray spectra, and an individual ULX can tran¬ 
sition between the two spectral modes (Kubota et al. 


2001), presumably due to variability in accretion onto 
the central source. Many of these objects in the local uni¬ 
verse (within 10 Mpc) are associated with lar ge ionized 
nebulae up to several hundred pc in diameter ( Pakull & 
Mirioni|2CI02| Cseh et al.|2012 Dopita et al. 2012). 'I'hese 
nebulae often show bright Ha bubbles, with the X-ray 
source either within the bubble or on one edge. Opti¬ 
cally, they are generally associated with strong narrow¬ 
line He 11 A4686A emission. This high-ionization line is 
not readily created by normal 0-stars, though broad He 
11 lines are observed in Wolf-Rayet stars. 



Figure 10. 3-color image of the central bubble. Background is 
the s ame as figurelTl right p anel. Chandra-identified X-ray sources, 
from Clements et m. ](l2002t. are mark ed as yello w cro sses, and the 
two iNlK nuclei l|fcicoville et al.||1998l see section |2.3|l are marked 
as white circles, i^tom east to west, the three sources shown are 
their X-4, X-1 (nuclear), and X-3 (X-2 is out of the field of view). 
The X-ray source coordinates have been shifted so as to align their 
X-1 nuclear source with the western NIR nucleus; thus the western 
nucleus marker is obscured. X-1 and X-4 are both unresolved hard 
X-ray peaks, while X-3 is the peak of extended soft X-ray emission 
(extent of 2”, shown as the yellow circle). 


The bubble feature in Arp 220 resembles these ULX 
bubbles. Its physical extent, 600 pc in diameter, is 
roughly the same size as the known ULX bubbles. It’s 
also coincident with a bright (Lo.s-io.ofeeV ^ 10®® ergs 




soft X -ray peak observed by Chandra (Clements 


et al. 2002, their X-3) (figure 10). However, while the 


X-ray sources associated with tne local ULX bubbles are 
generally compact, Arp 220-X-3 appears to be extended, 
or embedded in t he same extended emis sion that sur¬ 
rounds the nuclei. Clements et al. (2002) report an ex¬ 
tent of 2”, shown as the yellow circle in hgurefTOj No data 
exist on possible variability of this source, or of narrow¬ 
line He H emission. Finally, given the alignment of the 
bubble feature with the large-scale X-ray outflow, it is 
unlikely that the bubble feature has an origin unrelated 
to this structure. The large-scale outflow is bipolar and 
thus likely nuclear in origin. Further study is needed 
to explore whether this X-ray source is indeed a ULX 
with an affiliated nebula, or whether the X-ray emission 
is from hot gas filling the bubble created via one of the 
mechanisms described above. 


4 . 3 . Bubble origins, summarized 

We have explored four possibilities for the origin of the 
bubble: an outflow from an obscured AGN, an outflow 
from high levels of star formation within the nuclear gas 
disk, a traditional superbubble surrounding an undiscov¬ 
ered star cluster at the bubble’s center, or a ULX bubble. 
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Given the alignment of the nuclei, bubble feature, and 
large-scale (biconal) outflow structure, as well as energet¬ 
ics arguments, it is less likely that the feature is caused 
by an off-nuclear source unrelated to the outflow. While 
a chance alignment of a ULX superbubble is possible, 
and warrants further study, it is more likely that an out¬ 
flow or jet from the nuclear region is driving both the 
bubble and the outflow. 

Out of the two nuclear origin scenarios, we cannot dis¬ 
tinguish between an AGN jet or outflow or a starburst 
outflow. The low luminosity of the suspected AGN in 
the western nucleus (~10^^ erg s“^) and the high lev¬ 
els of star formation within the inner nuclear gas disks 
around both nuclei would each provide a similar rate of 
energy injection, either of which would be sufficient to 
drive the bubble feature. It is also possible that both 
mechanisms are contributing to the outflow driving the 
bubble. 


5. CONCLUSIONS 

We have resolved the near-nuclear Hq;-|-[NII] emission 
peak into a bubble feature using HST/WFG3 narrow- 
band photometry. Our new multi-wavelength imaging 
maps connect this feature to both the nuclear molecular 
gas disk and large-scale outflows. NIR IFS data show an 
increase in velocity dispersion at the base of the bubble, 
while the optical IFS data show that the emission in this 
region is dominated by shock excitation. This strongly 
indicates that the bubble is formed by an interaction of 
a jet or outflow, originating in or near the nucleus, with 
the surrounding ISM. This outflow shocks as it propa¬ 
gates into the surrounding galaxy and is connected spa¬ 
tially with the large-scale outflow that extends for multi¬ 
ple kpc from the galaxy. We discuss four possibilities of 
the origin of the bubble and find it most plausible that 
the bubble originates within the inner ~I00 pc around 
the nuclei. Either an obscured AGN or the observed high 
level of nuclear star formation could power the observed 
feature. Further work is needed to differentiate between 
these two possibilities. 
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APPENDIX 

Narrowband imaging of [SII]AA 6717, 673lA, in con¬ 
junction with similar im aging of Ha, has been used previ¬ 
ously on dwarf galaxies ( Galzetti et al.|2004 ) to separate 
photoionized from non-photoionized gas, and hence to 
evaluate photoionized gas oxygen abundances and non- 
photoionized gas distributions. One of the primary moti¬ 
vations for obtaining the Arp 220 HST dataset presented 
in this paper was to apply this same technique on a late- 
stage merger. Unfortunately, the [SII] emission line flux 
could not be recovered from the E680N filter, due to 


both complicating factors specific to the combination of 
the redshift of Arp 220 and this filter, and to the shal¬ 
lowness of the image. In this Appendix, we describe the 
reduction process for this filter and suggest ways to im¬ 
prove the process. 

The E680N filter, which contains the [SII] lines, is also 
contaminated by some of the Hq!-|-[NII] line complex, 
which falls on the blue cutoff of the F680N bandpass 
(figure 1^ . Because the bandpass cutoff falls in the mid¬ 
dle of the intruding line complex, we cannot simply sub¬ 
tract the flux-calibrated Ha -I- [Nil] image from the [SII] 
image. Instead, we must estimate the percentage of the 
contaminating lines that fall in the bandpass. We con¬ 
volved the SDSS spectrum with the theoretical through¬ 
puts of the F665N and the F680N filters obtained from 
the synphot PyRAfpd package to calculate the contri¬ 
bution of the flux of these lines. Roughly half the flux in 
the Had-[Nil] line complex in the SDSS spectrum falls 
in the E680N bandpass, though the exact contribution 
depends on the width and velocity of the lines, metallic- 
ity, and gas excitation mechanism. This will likely vary 
across the face of the galaxy, and introduces additional 
uncertainties. 

Ghanges in the line of sight velocity may shift the 
Had-[Nil] lines in or out of the F680N bandpass, altering 
the scaling factor used to remove the flux contributed by 
these lines. This error was estimated by red and blue 
shifting the SDSS spectrum by 500 km s“^ along the 
line of sight. At the velocity extremes, the flux contribu¬ 
tion from the Had-[Nil] lines complex ranges from 35% 
to 80% of the flux in these lines. The additional error 
in flux is included in the Icr errors on the [SII] image in 
table [TJ 

However, subtracting 50% of the Had-[Nil] flux from 
the continuum-subtracted [SII] image produced oversub¬ 
tracted images with negative flux values near the center 
of the galaxy, clear artifacts of the Had-[Nil] subtrac¬ 
tion. This is likely due to high extinction in the central 
regions. We explored several avenues for improving con¬ 
tinuum subtraction in this filter, all of which were ul¬ 
timately unsuccessful at producing robust detections of 
[SII]. This is likely due to the intrinsic faintness of the 
[SII] lines, which are only 10% of the flux in Had-[Nil] 
on average, and the short exposure time in the E680N 
filter. 

To make use of the [SII] data, we would require a longer 
exposure time in the E680N filter to allow us to probe 
the fainter average [SII] emission. Alternatively, separate 
high spatial resolution map of both [Nil] and Ha flux 
and velocity structure would allow us to more robustly 
correct the [SII] image for the [Nil] contamination, and 
hence recover more [SII] flux. Without either of these, 
the current [SII] data are inadequate for our purposes. 
Future projects that attempt this type of analysis should 
be especially careful to obtain deep images of the fainter 
emission lines. 
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